An Instability-driven Dynamo 
for Gamma-ray Burst Black Holes 



Rafael Angel Araya-Gochez 



Laboratorio de Investigaciones Astrofisicas, Escuela de Fisica 
Universidad de Costa Rica, San Jose, Costa Rica. 



Abstract. We show that an MHD-instabiUty driven dynamo (IDD) operating in a hot 
accretion disk is capable of generating energeticaUy adequate magnetic flux deposition 
rates above and below a mildly advective accretion disk structure. The dynamo is 
driven by the magnetorotational instability (MRl) of a toroidal field in a shear flow 
and is limited by the buoyancy of 'horizontal' flux and by reconnection in the turbulent 
medium. The efficiency of magnetic energy deposition is estimated to be comparable 
to the neutrino losses although an MHD coUimation mechanism may deem this process 
a more viable alternative to neutrino-burst-driven models of 7-ray bursts. 



I INTRODUCTION 

The combined redshift and fluence measurements of at least five 7-ray burst 
sources plus very large photon energy detections in certain bursts and tight size 
constraints derived from the rapid risetimes of burst triggers strongly suggest that 
the release of energy is highly focused by the central engine that propels a 7- 
ray burst. In spite of the very large 7-ray energy requirements, the efficiency of 
energy deposition into electromagnetic channels is likely to be very poor if the 
burst is driven by a neutrino burst in analogy to the processes thought to give 
birth to supernovae (MacFadyen & Woosley 1998, 1999) or if the burst involves 
major energy losses to gravitational radiation such as might be the case in compact 
object merger scenarios (Rasio and Shapiro 1994, Davies et al. 1994, Ruffert & 
Janka 1998). Indeed, these measurements pose a serious energy budget problem 
for arguably all gravitational collapse powered models of 7-ray bursts if the energy 
release is not moderately collimated. 

An attractive solution to this problem starts with a (directional) Poynting-flux 
dominated outflow (Thompson 1994, Meszaros & Rees 1997) under the premise 
that such a flow may carry very little baryonic contamination if deposited along a 
centrifugally (or gravitationally) evacuated funnel such as the angular momentum 
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axis of a black hole- accretion disk system. Yet, formal motivation for an external 
field of the desired strength and topology has yet to be investigated in this setting. 

Wc motivate a reasonable set of heuristic two dimensional dynamo equations for 
magnetic field components in the comoving frame of a mildly advective disk under 
the premise of negligible generation of meridional field (invoking Parker's undulate 
instability to promote the growth (loss) of vertical (horizontal) field is questionable 
in turbulent disks where the turbulence is fed by MRI's). A self-consistent turbulent 
steady state is achieved when the non-linear damping rate of the turbulent cascade 
equals the inverse of the linear growth timescale (Zhang, Diamond & Vishniac 
1994). Accretion disk models are far from being self-consistent in this respect. 

The heuristic dynamo equations account for field generation by the shear flow and 
by the non-axisymmetric MRI in Paczyhski's pseudo-potential well $ oc (r — Vg)^^ 
where Vg = 2GM/c^\ with flux buoyancy and turbulent reconnection providing 
for field loss terms. The model predicts azimuthally averaged field components 
and depends explicitly on the magnetic Mach number of the turbulence (assumed 
Alfvenic), on local pressure ratios, and on the relativistic generalization of the shear 
parameter (Aoort — Oort's A constant). 



II DYNAMO EQUATIONS AND SCALINGS 

In Lagrangian coordinates 

0,B^^^-^-^ and dtBr^^-^-^. (1) 

^5 Tm '^T 

The shear flow is parameterized linearly by a relativistic generalization of Oort's 
first constant (Novikov & Thorne 1974), = 2Aoort — T^c^inr^^, where 7 is the 
bulk Lorentz factor of the flow. Shear forces the radial wavenumber of perturbations 
to evolve according to kr{t) = — 2Ak^t. 

We use exact analytical MRI scalings from Foglizzo & Tagger (1995). With 
= B^/{87rp) and A = A/Q < 0, the maximum growth of non-axisymmetric 
MRI modes occurs for wavenumber 

K^l^ '^'{2A + (l+a,)A'}^ (2) 

at a rate, r]^^ that obeys 

|A|Vi"-^'l + «,(2 + A) (3) 

as long as k^^l/r. On the other hand, in a very strongly sheared flow 

-A > 1 + {2a^ + iy\ (4) 

the slow branch of MHD propagation (to which the toroidal MRI belongs) is desta- 
bilized into a radial interchange mode (T. Foglizzo, Priv. Comm.) in accordance 
with the Rayleigh criterium. 



A plausible mechanism that promotes baryon unloading from field lines is tur- 
bulent pumping (Vishniac 1995a) by the MRI which must favorably compete with 

turbulent diffusion of matter back onto flux ropes. Under this assumption, the 
stretch, twist and fold of field lines by (enthalpy- weighted) sub-Alfvenic turbulence 
augments the field energy density and releases matter from field lines that other- 
wise would be "frozen-in". For the marginal case of Alfvenic turbulence, nearly 
empty flux ropes (i.e. flux residing on surfaces perpendicular to the local 
meridian) in a gas pressure dominated disk acquire a drag limited buoyant velocity 
Vb oc {v^fY /Cs- Moreover, assuming efficient diffusion of radiation and pairs into 
the flux tubes, in this picture the buoyancy loss rate, ^ = Vh/He, is enhanced by 
a factor s p/pgas = C (Vishniac 1995b). 

Reconnection of the field at sub-MRI optimal lengthscales (where the field lines 
are only weakly stochastic) occurs at the Alfven speed for Alfvenic turbulence 
(Lazarian & Vishniac 1999). The rates may be written as inverse Alfven transit 
times calculated from the component of the field that undergoes reconnection. 

Tree ~ l^Jvt^^ = yjojv {I^I/Hq) (cs/w/^'O whcrc r is the adiabatic index of 
the fluid. 

The perpendicular lengthscales, associated with the mean distance for field 
reversal are derived from the fundamental linear lengthscales for coherent field 
pumping. These arc supplied by half of the toroidal MRl's (wave)length scale. 
Azimuthal, radial field reversal, 1^, directly involves the optimal wavenumber of the 
non-axisymmetric MRI. In addition, following Tout & Pringlc (1992), an estimate 
of the radial, azimuthal field reversal lengthscale, follows by noting that the 
time evolution of wavenumbers implied by shear during one MRI timescale couples 
the azimuthal lengthscale to the radial lengthscale, i.e. ly = (Tj^/t^) x l^. Thus 
Yr = tt/^m and = (t^/t^) x (tt/A;^). 



Ill EQUILIBRIUM SOLUTIONS 

Scaling wavenumbers to the inverse pressure scale height k — > k/He, and re- 
defining A = |A|/f2 > 0; a set of normalized dynamo equations follows by rep- 
resenting fields in (Alfven ) velocity units and normalized to the soundspeed 
(B' = B X y^AWgCs), and the time normalized to the inverse of the Keplerian 
frequency t' — D,t. 

In a steady state, these equations must satisfy 



which we solve numerically. 



IV THE ENERGY DEPOSITION RATE 



The accretion disk setting is envisioned to follow the standard hyper-accreting 
black hole model of Popham, Woosley & Fryer (1999, hereafter PWF) where ~ 
3 M©; ttss = 0.1, and M = 0.1 sec^^ By adopting their published pressure 
ratios and assuming Keplerian rotation for r e [2.25, 20]rg, we "piggyback" the 
hydromagnetic energy conversion process on this model. 

We find that the magnetic output rate from buoyancy, E ~ 1.77+51 erg scc~^, 
is comparable to the neutrino luminosity Ly ^ 3.3_|_5i erg sec~^ (PWF). The 'half- 
luminosity' radius is located at ~ 5.75rg and the IDD becomes operational 
(against the radial interchange instability, c.f. Eq [4] and Araya-Gochez 1999) at 
rjnin — 2.55 r^. Curiously, the derived value of the magnetic viscosity ass = B'^B'^ 
hovers on 0.1 at the innermost radii and decreases to about .087 at rc in good 
agreement with the adopted value. Lastly, we note that the MRI pumps the field 
preferably at the lowest wavenumbers for a near-equipartition field in a strongly 
sheared flow. Thus, most of this energy could in principle go into a Poynting jet 
with interesting consequences for 7-ray bursts. 
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